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Abstract. The study of factoring relations between subshifts or cellular 
automata is central in symbolic dynamics. Besides, a notion of intrinsic 
universality for cellular automata based on an operation of rescaling is 
receiving more and more attention in the literature. In this paper, we 
propose to study the factoring relation up to rescalings, and ask for the 
existence of universal objects for that simulation relation. 
In classical simulations of a system S by a system T, the simulation takes 
place on a specific subset of configurations of T depending on S (this is 
the case for intrinsic universality). Our setting, however, asks for every 
configurations of T to have a meaningful interpretation in S. Despite 
this strong requirement, we show that there exists a cellular automaton 
able to simulate any other in a large class containing arbitrarily complex 
ones. We also consider the case of subshifts and, using arguments from 
recursion theory, we give negative results about the existence of universal 
objects in some classes. 



1 Introduction and definitions 

Tilings and cellular automata are two paradigmatic models often considered 
in the fields of complex systems and natural computing. They are comple- 
mentary — one is static and non-deterministic and the other is dynamic and 
deterministic — but they are both formally simple and both related to symbolic 
spaces. Moreover, many links are now established between the two models (see 
for instance [10 9 ) so it is natural to consider them together. 

Both are known to be Turing-powerful since their introduction in the mid- 
20th century [24118) . However, analyzing their ability to process information 
only through translations into the Turing world is very restrictive. Such models of 
natural computing deserve a natural and intrinsic notion of reduction to compare 
their objects one to each other. Following this line of thought, several notions 
of simulations were proposed recently which are intrinsic to each model, and 
lead to corresponding intrinsic notions of universality |19I23I15T5| : a system is 
universal if it is able to simulate any other from the same class. 

Intrinsic universality for cellular automata is probably the most studied of 
such notions |20|21|16|2|6j . The underlying relation of simulation uses uniform 
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encodings working at the level of blocks of cells. More precisely, S simulates T if 
S, when restricted to a suitable subset of 'correct' configurations, is isomorphic 
to T via such an encoding. Our approach is different and uses redundancy of 
information instead of restriction to a subset of configurations. In our setting, 
S simulates T if there is a uniform way of projecting the whole phase space of 
S onto the phase space of T (a precise definition is given below). The question 
addressed by this paper is the existence of universal objects with respect to that 
simulation relation, we call them factor universal objects. 

The first contribution of this paper is the formalism based on the well-known 
mathematical notion of action: it allows to encompass both subshifts and cel- 
lular automata, it gives a new look at the notion of cell grouping which is the 
root of the simulation relation used in intrinsic universality, and it establishes 
connections with the work of Hochman [S] where the use of sub-actions is crucial. 
Our main result is that, although factor-universal objects do not generally exist 
(theorem [IJ, it can still be constructed for some large class like the set of cellular 
automata having a persistent state (theorem [2]) . 

Basic definitions. Given a finite set Q and an integer d > 1, the symbolic space 
of dimension d over alphabet Q is the set . It can be seen as an infinite set 
of cells arranged as a lattice Z'' and each carrying a value from Q. An element 
of is called a configuration. is naturally equipped with the compact 
Cantor topology [T3] which is the product topology of the discrete topology on 
Q (it can also be defined via a metric). 

Another key notion in the context of symbolic spaces is that of finite patterns 
that may occur in infinite configurations. For our purpose, rectangular patterns 
will be enough. Given z = {zi, ■ ■ ■ , Zd) G Z'^ with > for all i, the hyperrect- 
angle TZz is the set of vectors z' = {z{, - ■ ■ , z^) G such that Q < z[ < zi for 
all i. A Q-pattern of shape TZz is a coloring of TZz by Q, that is an element of 
Q^'' . Given a configuration c G , the pattern of shape TZz^ extracted from c 
at position Zp e lA, denoted by T'f^ (c), is simply: z e TZz^ >-> c{zp + z). 

The objects we study (subshifts and cellular automata) share the property of 
being uniform, i.e. invariant by translations. Formally, given z G Z'' the trans- 
lation of vector z, denoted Uz, is the function mapping a configuration c G 
to the configuration az{c) such that Vz' G Z'', az{c){z') = c{z' + z). 

A subshift is a subset of which is translation invariant and closed for the 
Cantor topology. Equivalently, a subshift is a set El of configurations avoiding 
any occurrence of any finite pattern from a given language of patterns L: 

= {c G g^" : Vz, z' G Z'' with z, > for aU i, V^, (c) ^ L). 

A subshift of finite type is a subshift of the form where L is finite. There 
are strong connections between subshifts of finite type in dimension 2 and sets 
of tilings generated by a set of wang tiles. In particular, due to Berger's theorem 
[I], it is undecidable, given a finite L, to determine whether E^is empty or not. 

A cellular automaton is a local and uniform map on a symbolic space. For- 
mally, it is given as a 4-tuple by its dimension d, its alphabet Q, its neighborhood 



1/ C (finite) and its local transition map f : — s> Q. To that formal object 
we associate a global map F acting on as follows: 

Vc e Z'', Vz e Z^ i^(c)(z) = /(z' ^c{z + z')) . 

The fundamental theorem of Curtis-Lyndon- Hedlund [7] states that global maps 
of cellular automata are exactly continuous maps on symbolic spaces which com- 
mute with translations. 

Actions and rescalings. Let (M, +) be a monoid (a set equipped with an 
associative law and a neutral element). An Ml- action on a space X is a function 
\I/ -.Mx X ^ X such that !f'(0, x) = x (for &\\ x € X and being the neutral 
element of M) and 

Vx e X, Vm, m' G M, !f'(TO + m', a;) = S'(to, ^{m' , x)) . 

We will use the formalism of action to study both subshifts and cellular au- 
tomata: 

— if C is a subshift, we canonically associate to it the Z''-action on 
U defined by Ws{z,x) = az{x); 

— if is a cellular automaton on the space , we canonically associate to it 
the N X Z''-action on Q^" defined by <I'F{{t, z),x) = cr^ o F\x). 

If M' is a sub-monoid of M, 'F induces an M'-action by restriction to the 
domain M' x X. M and M' can be isomorphic or not and both cases might be 
interesting. For instance, studying a cellular automaton _F as a classical dynam- 
ical system consists in forgetting the spacial component of Fp and focusing on 
the pure temporal action of F. This point of view was often adopted in the lit- 
erature (e.^., topological dynamics of cellular automata [13 ) but, interestingly 
enough, recent work of Sablik [22] tends to re-incorporate the spacial component 
of actions to better study the dynamics of cellular automata. 

In this paper, we will only consider the case where M and M' are isomorphic. 
More precisely, in our context, M will be of the form Z*^ or N x Z'' and we 
will consider sub-monoids of the form M' = toN x ziZ x • • • x Zd'Z, with tg > 
and Zi > for all i. In this case, passing from the M-action to the M'-action 
can be seen as a neutral change of point of view on the system that we call 
rescaling in the sequel. The intuition is that we change the discrete units of 
time and space, passing from 1 to to in time and 1 to Zi in direction i. Given 
a subshift or a cellular automaton, a scaled action is simply the restriction of 
their canonical action to some sub-monoid of the form M'. It is worth noticing 
that a scaled action associated to a subshift (resp. a cellular automaton) on 
the alphabet Q is always isomorphic to the canonical action of a subshift (resp. 
a cellular automaton) on an alphabet of the form . More concretely, this 

isomorphism comes from the natural one-to-one map from to (Q^^")^ , 
where Zs = (zi, . . . , Zd), which maps a configuration c to: z i— )> V^'^^ (c), where 
the operation x on 1/'' denotes coordinate-wise multiplication. Our notion of 
rescaling for cellular automata is similar to the one in [19123) which is the basic 
ingredient to define intrinsic universality. 



Factors. One of the central notion in symbolic dynamics is that of factor. 
Intuitively, a factor is a uniform continuous projection. This notion has also 
been used with success in the study of expansive cellular automata |17j and more 
generally as a classification tools for cellular automata |12I3| . As we study both 
multi-dimensional subshifts and cellular automata, we give a unified definition 
using the formalism of actions. 

Definition 1 Let M and M' be isomorphic monoids via i : M M'. We say an 

Wl' -action (j)' on X' is a factor of a W.- action ip on X if there is a continuous onto 
map TT : X ^ X' such that: Vx G X,\fm G M, Tr(if){m,x)) — (f)' (i{m),Tr{x)) . 

Two key points are that: (1) any orbit in {(l),X) projects onto some orbit of 
{(j)',X') via TT, and (2) any orbit of (p' can be realized as such a projection. 
In a word, the simulation of {(/)', X') by {(l),X) is everywhere meaningful and 
complete. 

2 Factor Universality 

At this point, we could compare subshifts or cellular automata through the fac- 
toring relation between their canonical actions, saying that system S factors 
onto system T if the canonical action of S factors onto that of T. However, this 
gives an excessive importance to the alphabet and forbid the existence of uni- 
versal objects due to entropy considerations (factoring cannot increase entropy) . 
In f5|, this limitation is bypassed via dimension changes: a d-dimensional system 
is compared to fc-dimensional systems (fc < d) via its fc-dimensional sub-actions. 
Our point of view is different. We always work at constant dimension, but we use 
another kind of sub-actions: scaled actions defined above. For a fixed dimension, 
monoids of scaled actions are all isomorphic and we will consider only canonical 
component-wise isomorphisms between them. We can now formulate the central 
definition of the paper. 

Definition 2 Let S and T be two d-dimensional subshifts (resp. CA). We say 
that T is simulated by S, denoted T ^ S , if some scaled action of S factors onto 
some scaled action of T . 

As usual when working on symbolic spaces, continuity and uniformity implies 
locality (Curtis-Lyndon-Hedlund theorem [7]). In our context of rescalings, the 
locality is no longer expressed at the level of cells, but at the level of groups of 
cells. More precisely, we say that a map (j) : Qf — >■ Qf is local if there exist: 
r G N (locality radius) , two shapes TZ^^ and TZz2 (source and destination scales) , 
and a local function / : Q^^^'"''^-'^'^ — s> Q^^^ such that 

vc G gr,vz G (0(c)) ^ f{v^^;t:t\^ (c)). 

To fix ideas, ii d — zi — Z2 — 1 and Qi = Q2, f just the local map of a 
cellular automaton of radius r and (j) is its corresponding global map. 



Proposition 1 Fix a dimension d. Let Si and S2 be two d-dimensional sub- 
shifts and let Fi and F2 be two d-dimensional CA of alphabet Qi and Q2 respec- 
tively. Then we have: 

— E2 ^ Si if and only if there is a local map such that 4'{Si) = S2; 

— F2 =4 Fi if O'Tid only if there is an onto local map cj) from to and 
integers ii, ^2 £ N such that 4> o F*^ = F^^ o 0. 

Besides the work of Hochman [8 , notions of simulations similar to ^ have 
already been considered for tilings jl5j or for cellular automata [2314] or more 
general settings [13]. Each time, one of the main concern is the existence of 
universal objects: this is precisely the central point of the present paper. 

Definition 3 Let C be a class of sub shifts (resp. cellular automata). A subshift 
(resp. cellular automaton) U is C-universal if U ^ C and X ^ U for any X eC. 

Whatever the fixed dimension, there is no universal subshift for cardinality 
reasons: there are uncountably many subshifts but for a given subshift U there 
are at most countably many different subshifts =^-simulated by U (by proposi- 
tion [1]). The following theorem uses recursion theoretic arguments to yield other 
negative results concerning universality (similar arguments where used in |19I8) 
in different settings). 

Theorem 1 Fix a dimension d > 2. Then there is no universal subshift of finite 
type of dimension d and there is no surjective- universal CA of dimension d. 

3 A Large Class with a Universal Object 

In this section, we restrict to dimension 1 to make a clear exposition of the main 
result (theorem [5]) . 

Definition 4 A CA A is said to be persistent if there is a state Qq £ such 
that for any configuration c £ Qa" if = go then A{c) (i) = go- 
We denote by P the set of all persistent CA. 

Note that for any CA, you may add an extra persistent state and obtain a 
CA in P containing the dynamics of the first one. 

Theorem 2 There exists a P-universal cellular automaton. 

In the following, we describe such a P-universal CA U with radius 1 and 
alphabet Qu. 

We denote by Pq the set of CA of P with radius 1 and alphabet of size 2^ 
for some p E N. One may easily describe for any CA B € F a CA A € Fg 
such that B =4 A. Using transitivity of ^, it will be suflficient, in order to prove 
P-universality oi U, to exhibit for any CA ^ e Pq an onto local map (j>j, from 
Qu^ to (5_4^ and an integer ta such that U^-* o = 0^ o ^. 

To do so, for each A, we introduce an integer 1^ and a dichotomy on words 
of Qu'^. 



— on the one side we have what we caU A- correct macrocells (or ^-macrocells). 
They encode information about a current state x G Qa, about the local rule 
of and a machinery used to apply this rule to update the current state. 
In almost any case, they will be interpreted through (p^. as x. 

— on the other side we call all the other patterns A-incorrect, and they will be 
interpreted as the persistent state of A. 

The idea behind the local rule of U is to make every ^-macrocell determine 
if it is surrounded by other ^-macrocells. If this is the case, then interaction 
is possible, and the current state of the neighbor will be taken into account 
to compute the new current state, following the rule of A. Otherwise, there is 
no interaction, the .A-macrocell evolves considering every ^-incorrect neighbor- 
hood as a persistent state neighbor. The difficulty is that although correctness 
is related to the particular CA being simulated, every configuration must evolve 
correctly for every possible CA. 

The proof of universality uses the combination of two key properties: on one 
hand, correct patterns remain correct and evolve according to the rule being 
simulated, even if not surrounded by correct patterns (lemma S]); on the other 
hand, incorrect patterns are interpreted as the persistent state and never become 
correct (lemma [S|). 

To make the construction of lA readable, we describe its state set as a su- 
perposition of several layers: the main layer Ad contains most of the information 
about the simulation and the macrocells informations; signals layers are used to 
manage the evolution of the main layer; and clock layers guarantee synchroniza- 
tions. 

Correct macrocells description In the following we consider a simulated 
CA ^ e Pq with radius 1 and state set of size n — 2^. We use a canonical 
binary enumeration of the state set, in which the first word (0'°^^"^) represents 
a persistent state of A, denoted pA- Our ^-correct macrocells will be words of 
length I A (specified above) whose main layer follow the pattern 

# Ci \Transition table\ \State\ \memory\ # 

— # are delimiters, they never appear or disappear during the computation 

— Ci is the control state used to control the successive steps of computation. 

— \Transition table\ is the binary description of the transition table of A. 

— \State\ contains the binary value of the current state of the macrocell. 

— \memory\ is a binary area which will be used to keep the values of the neigh- 
bors' current states before computing the new current state of the macrocell. 

\Transition table\, \State\ and \memary\ are encoded with disjoint binary 
sub- alphabets. A cell whose state belongs to one of those sub-alphabets will 
never change sub-alphabet. Moreover, the states of the transition table's cells 
are never modified. 

The current state description is Zo(7(n)-bits-long. In the transition table, im- 
ages are ordered canonically, so the length of the description is simply n^log{n). 



The memory should be at least 2/05 (n)-bits- long in order to contain current 
state values of the two neighbors. But in order to simplify some proofs, we chose 

such that it is at least half the total size of the macrocell, and such that the 
function A ^ Ia one-to-one. 

Most of the computation will happen on those very constrained patterns. In 
the next definition, we add an extra constraint on the control state to obtain 
^-correct macrocells. 

By stability of the sub-alphabets, such a correct ^-macrocell remains correct, 
but the state value may change. This is why we take some care when defining 
the current state value associated to the macrocell. 

Definition 5 A word u € Qu^ of length Iji is said to be a ^-correct macrocells, 
denoted by u G C^, if its main layer follows the structure defined above ( correct 
sub-alphabets for each cell, and correct transition table of A), and if its control 
state is in Co- 

For each such A-correct macrocell u, we define its associated state value 

v{u) £ (3_4, which is the state described by its current state value after steps 
of computation by U. And this state value v{u) only depends on u. 

To ensure that the sate value only depends on w, the memory area is used 
as a buffer to prevent modification of the current state value coming from the 
left, before the computation has been initialized. Details will be given in the 
following. 

By extension we may sometime call ^-macrocells words following the general 
pattern, even with non-Co control state, in particular when they are images of 
a ^-correct macrocell. 

The local rule. We describe the local behavior of U starting from a correct 
vA-macrocell. The local rule will first determine which neighbors it may interact 

with {Check of the neighbor's length and synchronization, and Transition table 
and state encoding check), and then compute its new current state according 
to the rule of A and eventually the value of those neighbors {New current state 
computation). 

In order to guarantee the synchronization between ^-macrocells, we specify 
the duration of each step, and even of some sub-steps. It is done by a clock, 
which use specific layers of the states; their existence is proved by the following 
lemma: 

Lemma 1 For any k,h £N \ {!}, there exist a CA, and two states q^, and qf 
such that the leftmost cell of an area delimited by two # separated by I — 2 cells 
turns to state qf at some time t > k-P -\-h-l iff this cell was in state qs exactly 
k ■ P + h - 1 steps before. Moreover, this property is guaranteed independently of 
what is outside the two #. 

At the beginning of each step, the control state C, will turn to Cj+i, ini- 
tiate the corresponding clock, and initiate some signals which will manage the 



evolution. Those signals are distinct states propagating on upper layers of the 
configuration, and interacting with the main layer and other signals. We say 
that a signal belongs to a macrocell if it was generated in this macrocell's area, 
between the two And, thanks to our evolution rule, a signal always knows 
if it is in its cell or in the area to the right or left of its cell. It is also useful 
sometimes to make signals carry one extra bit of information. It is simple to do 
it using distinct states, since the number of bits is bounded. 

Check of the neighbor's length and synchronization. Co — > Ci (that is 
to say that when the control cell's state is Co it becomes Ci): 

Recall that we are interested to the behavior in the case of an ^-correct 
macrocell. When Cq becomes Ci, it initializes two control bits with value 0, in 
the main layer of the control cell, and it launches signals. Since the construction 
is classical, we simply illustrate the desired behaviors by figure [TJ Those two 
pictures illustrate the signal machinery in the case of respectively left and right 
neighbors of same length and with state Cq appearing simultaneously (what we 
call synchronized) . Every transition whose image is one of the signal involved in 
this checking appears on those pictures. 

The first signals Si and S4 erase all signals belonging to our macrocell. To- 
gether with the length of the memory being bigger than the length between the 
control state where signals are generated and the end of the current state area, 
it constitutes the protection of the current state value from eratic signals coming 
from outside the macrocell, and justifies that in definition [5] w(u) is a function 
of u. 

If the neighbors have same length and are synchronized, this whole step takes 
4 times the length of the macrocell, 1^. After 4 • 1_a steps, the control state Ci 
becomes C2, and if it did not receive a positive result from one side, it concludes 
that the involved neighbor is incorrect. This is managed using a clock signal 
(with h = 2 and fc = in lemma [Ij initialized by Cq on a specific layer. When 
Qf is raised on this layer, Ci becomes C2. The important point is that we ensure 
the following property. 

Lemma 2 The control state of an A- correct macrocell becomes C2 exactly 4 • Z_4 
steps after Co appeared. At this step each control bit has turned to 1 iff the 
corresponding neighboring macrocell has same length and synchronisation than 
the considered macrocell. 

The proof of this lemma is direct for the length but asks to enter into some more 
(simple but fastidious) details for the synchronization part. 

Transition table and state encoding check. C2 — !> C3 : 

In this step, for each neighboring pattern with same length and synchroniza- 
tion, the macrocell checks whether the transition table and the current state 
are compatible with its own (same lengths, and same content for the transition 
table) or not. 

When C2 appears it launches the following test for each neighbor whose 
corresponding control bit was 1, and initializes two fresh bits to 0. First, a signal 



Fig. 1. Successful left (resp. right) neighbor test by the right (resp. left) macro- 
cell (mix of main and signal layers for easier reading) 

is generated and puts a mark (that is to say a non-moving signal) on the first cell 
of the transition table of its macrocell, and another mark on the first cell of the 
transition table of the neighbor it checks. Then signals are exchanged between 
those two marks that will each time carry the binary state of the cell pointed 
by one mark to the next unchecked cell of the other macrocell; it compares this 
cell's binary state to the carried binary state, and push the mark by two cells. If 
no difference is detected and if both marks reach the end of the transition tables 
simultaneously, a correctness signal is sent to the control state. 

After the transition table has been checked, the same mechanism is used to 
check that the current state encoding areas have same length. At the end of those 
tests, the results are sent to the control cell which again keeps the information 
on two control bits. For each cell of the transition table or the current state, 
checking takes 2 • Ij^ steps. So checking a whole neighbor takes less than 2 • l^. 
Again, a clock is used to make this test last exactly 2 • steps. Then the control 
cell is turned to C4. 

Lemma 3 The control state of our macrocell becomes C4 exactly 2-1^ steps after 
C'2 appeared. At this step each control bit has turned to 1 iff the corresponding 
neighboring pattern has same length, synchronization, and if the lengths of the 
transition tables and current states, and the content of the transition table are 
equal. In this case we say that this pattern is compatible with our A-macrocell. 



The proof of this lemma is straightforward. Keep in mind that some signals 
erased all erratic signals that could interact with our cell at a previous step. 

New current state computation. C4 ^ Ci^ : 

After all the tests have been done, the new current state has to be computed. 
We need to explicit how we consider the neighboring pattern. In the following, 
what we call detected state of one such pattern by our macrocell will be: either 
the persistent state if the neighbor is non-compatible with our y^-macrocell, or 
its current state if this is a compatible ^-macrocell. 

At first, the detected states of the left and right neighbors are written to the 
memory. It is written in the binary memory alphabet. Each detected state is 
written on login) cells. If one neighbor is compatible, we copy its current state 
value to the memory using marks and signals similarly to the previous step. If 
it is not compatible, we write 0'°^'"-', the length being the same as that of the 
current state area. We add a clock to specify that copies last exactly 2 • steps, 
the neighbor being correct or not. 

After 2 • additional steps, the search for the image in the transition table 
starts. It consists in reading the binary word formed by the three image states 
(the current state of the cell followed by the two detected states copied in the 
memory), and turning it into a unary position in the transition table. It is then 
possible to place a mark at this position, and finally copy this pointed state to 
the current state area. We make the reading of the position last 4 • l\ steps. And 
copying the new state lasts 2 • l\. After the whole computation step, which lasts 
7 • /3^, the control state turns to C5. 

Finally one step of simulation is completed after exactly tj\^ = 9 ■ l\ + A ■ Ij^ 
steps. After this time the control state turns to C5. 

To become Cq again, and launch a new step of computation, we add another 
condition. We ask a clock launched exactly steps before to raise a fiag. And 
obviously this clock may only be launched by Cq. It is realized using again signals 
of the lemma [T] computing on one more layer. 

The state set of the universal CA is given by Qu = Af x S* x C U {C/} with 

— the main layer : M = {Co, C5} U {Ci}ig{i,..,4} x {0, 1}^ U {Oi, li}ie{tt,cs,m} 

— the signals layers : S = Xig/{s,;} x Xjgj({sj} x {0,1}) 

— the clocks layers (see lemma[T]), one for each duration needed. C ~ ({0, 1} x 

{Sj}»G/c)** 

— is a single persistent state ensuring that W G P 

Yet, the transition rule of U is partially specified, we call correct transitions 
those defined up to now, in the case of correct macrocells. But the other transi- 
tions may not be chosen arbitrarily. We specify the following behaviors: 

— C/ is never modified by any transition 

— the main layer is never modified by a non-correct transition: they act as the 
identity on the main layer. 

— concerning the signal layer, apart from the collisions corresponding to the 
behavior described in the previous steps, all signals may cross each other 



(each kind of signal is evolving on its own layer). However, except for tran- 
sitions involved in the behavior described above, any signal that crosses a # 
is destroyed. 

Interpretation We now describe the continuous onto map (pA '■ Qi/' ~^ Qa' 
associated to A. This map is induced by a local map from patterns of shape 
I A to individual states of A. More precisely, using notation from proposition [1] 
we have r = 0, zi = 1^, Z2 — ^, t2 — I and ti — t^- 

If PA is the persistent state of A, the local map 4'A is defined as follows: 

1. Vm ^ Ca, fA{u) = PA 

2. Vm S Ca, fAiu) = v{u), with v{u) the value from definition[5] 

Proof of theorem [2] The proof of the theorem relies on the two following 
lemmas. They are consequences of the construction, the intermediate lemmas 
and the clock lemma. 

Lemma 4 V c e Qu^, to G N, (c)[o,;^-i] G Ca, thenv = W*°+"^(c)[o,i^-i] e 
and tpA{v) = Sa{iIja{c[-i^.-i]),i1^a{c[o.i^-^i]),iI'a{c[i_^.2-u-i]))- 

Lemma 5 If 3 t > ta, c £ Qu^ such that u — Z//*(c)[o,i^-i] £ Ca then v = 
^*""^(c)[o,i^~i] eC^- 

We can finally prove our main claim: €fo,A^t{. We use the characteri- 
zation of proposition [TJ Let .4 S Po the associated length I a and function (j>A are 
defined as explained before. First, (j)A is local (by definition) and onto, because 
correct macrocells are enough to encode any state of A and thus concatenations 
of correct macrocells allow to encode any configuration of A. Second, we have 
(j)A°l^^^ — Ao(j)A. To see this we discuss on the pattern of shape TZi^ at position 
and the rest follows by translation. If this pattern is not in Ca its image after 
TA steps remains out of Ca (lemma [5]) . If conversely this central word belongs 
to Ca, lemma m gives the desired property. 

4 Perspectives 

A natural extension of our work could be to generalize the construction to cellular 
automata having an equicontinuous point. The idea would be to use blocking 
words as a replacement for the persistent state. But it seems much harder, if not 
impossible. 

Besides, the main open question left by this paper is the existence of universal 
CA. We conjecture that they do not exist and more precisely that no CA can 
simulate all products of shifts. A possible way to obtain this negative result would 
be to study limit sets: by a compacity argument, one can show that a universal 
CA must have a universal limit set. The main obstacle is that subshifts that are 
limit sets of CA are not well characterized. 

Finally, we also leave open the existence of universal SFT and universal 
surjective CA in dimension 1. 
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A Proofs from section [2] 



Proof (Proposition\^ . First, an onto local map from Si to S2 with shapes TZz^ 
and TZz^ induces a factoring relation from the M-scaled action of Si onto the 
M'-scaled action of S2 with 

M= (zi)iZ X ■••(2i)dZ 

and 

M' = (22)1^ X •■•(22)dZ. 

Conversely, suppose that the relation Si =<; S2 is realized by a factor map tt 
from the M-scaled action of Si onto the M'-scaled action of S2 with 

M= (2i)iZ X 

and 

M' = (22)1^ X •■ .(22)^2. 
Consider now each pattern j> G Q2 • Since the cyhnder Cp defined by 

= {c e Qf : Po^^ {c)=p} 

is both open and closed, so is 7r~^(Cp). By compacity, and since cylinders form 
a basis of the topology, we get that 7r~^(Cp) is a finite union of cylinders of . 
We can suppose without loss of generality that they are all of shape Ti(2r+i)zi 
for some large enough r (finite unions of cylinders of small shape can always be 
defined as finite unions of cylinders of larger shapes) . Doing this with the same 
value of r for all p, we get a (possibly partial) function / from (^^-H-i)"^! 
Q2 ^ ■ By eventually completing / and by definition of the factoring tt between 
M and M'-scaled actions, / induces a local map from Qf to Qf associated with 
shapes TZz^ and TZz^ - It is onto because tt is onto. 

For cellular automata, the reasoning is similar and adding the temporal com- 
ponent in actions translates exactly into the desired property of weak commu- 
tation between the global maps of cellular automata and the onto map between 
configuration spaces. □ 

Proof ( Theorem Ul). For the case of surjective CA, it is enough to notice that 
surjectivity is preserved by the relation Indeed, if F =^ G we have 

o G*i = F*^ o 4, 

for some onto map 0. Therefore F must be surjective if G is surjective. 

Then the proof follows from Kari's theorem [TT] establishing that surjective 
CA are not recursively enumerable. Indeed, given a surjective universal CA U, 
we could enumerate thanks to the local presentation of factors (proposition [T|) 
all CA F such that F =4 U: they are all surjective (surjectivity is preserved by 



factor) and all surjective CA are aniong them (universality). 



We consider now the case of subshifts of finite type. Without loss of generality, 
any subshift of finite type can be presented as a subshift El where i is a finite 
set of patterns having all the same shape TZz for some z. By Berger's theorem 
[I] the set of such L verifying that Sl is not empty can not be recursively 
enumerated. We show below that the existence of a universal subshift of finite 
type implies the existence of an algorithm of enumeration of all L of the form 
above such that Sl is not empty. 

So suppose that there exists some universal subshift of finite type Elu where 
Lu is a set of Qy-patterns of shape TZzu ■ Obviously, must be non-empty. 
For any L and any pattern p of larger shape, we say that p is i-valid if it 
contains no occurrence of any pattern from L (occurrence requires that one 
shape is completely included into the other). 

Let L he a set of Q-patterns of shape TZz and -0 be a local map from Qf^ 
to Q"^ associated to shapes TZz-^ and TZz^- Consider the minimal shape TZz_y_ 
containing both TZz and 7?.2z2 • Since ip is local, one can check in finite time the 
following property called validity property: any pattern p of shape TZzj^ which 
has a L(7-valid preimage via if) is L- valid (the size of preimages of finite patterns 
depends on the radius r associated to V but details don't matter here). By the 
definition of local maps and the hypothesis on shapes, this property implies 
that V'(-£'L[/) — and therefore Sl ^ % (the choice of shape 7?.2z2 ensures 
that validity is checked inside blocks of shape TZz^ but also across the boundary 
between two such adjacent blocks). 

It follows that we can recursively enumerate couples {L,ip) having the prop- 
erty above. More precisely, maps tp are enumerated via their local presentation 
(shapes, radius and local function). This way, we can enumerate a list of finite 
languages L such that El is not empty. To conclude the proof it is sufficient to 
show that all L such that El ^% are present in the list. Suppose by contradic- 
tion that some L over alphabet Q with El 7^ is such that no local map from 
to Q'^ verifies the validity property above. By universality of Elu, there 
exists a local map ip sending El^ to El- Let r, TZz-i, and TZz^i, and local function 
/, be the parameters associated to ^. For any A; > we can define the same 
map ij) with another presentation by increasing artificially the radius r to kr 
and changing the local function / accordingly (shapes are kept unchanged). We 
call it the k*'' presentation of ij}. Since, by hypothesis on L, no such presentation 
has the validity property, we deduce that there must exist some finite pattern p 
which is not L- valid and such that, for any fc, p has a Lt/-valid preimage under 
the k*'' presentation of i). Therefore, by a simple compacity argument, there 
exists c S Elu such that ^-!{c) has an occurrence of p. Hence, V'(c) ^ El which 
is a contradiction. □ 



B Proofs from section [H] 



Proof of lemma[Ji We first build a CA that satisfies our lemma for k — 2, h = 0. 
Its state set will be made of one binary layer, and a signal layer. The behavior is 
simple: when Qs appears it generates a signal that will keep oscillating between 
the When the signal is generated for the first time, it initialize the area, 
turning the first binary cell to 1 and the other one to Os. Then, the signal keep 
moving from right to left and back between the #. Each time it goes to the right, 
it turns one more binary cell to 1. And when the rightmost cell's binary layer is 
finally turned to 1 a new special signal is sent to the left which will generate the 

If two or more signals crosses, one of them may survive. If one of them is 
initializing it will survive. 

So, in 2.1 steps of computation, the total number of Is may be non increasing 
only in the following cases: 

— if there is no signal at all 

— if all cell's binary layer is already 1 and in this case a qf was generated 

— if an initialization signal has been sent. 

In particular, if a qj appears at some step, then in the previous 4.Z steps, a 
1 was generated. And in each previous 2.1 step, at least a 1 was generated. 

Thus, in the previous 2.P steps, at least one initialization signal was launched 
and a qs has appeared. But by construction, after a qs state appears, the first 
qf state appears only exactly 2.1"^ steps later. 

It concludes the proof of the clock lemma in case k = 2, h = 0. For other 
values of k, simply slow down the signal going right to left. For other values of 
h, after the end of the quadratic part, launch a signal that will go right with 
speed 1 and come back left with speed — 1) before raising qf. □ 

Proof of lemma [5l' 

By definition of a correct pattern, u is given by: 

# Co \Transition table\ \State\ \memory\ # 

First of all, the # are never created or destroyed. The transition table and 
maximal state information are never modified, so they are the same in u and 
in V. And the sub-alphabet corresponding to control state, current state and 
memory alphabets are stable. The structure of v is the same as this of u. To 
prove our lemma it remains to prove that the second letter in w is Co, and that 
the current state value is smaller than the maximal value. 

But, to make Co appear, at step t, a signal qf was raised by the global clock, 
which implied, using the clock lemma, that the second letter of v is Cq. 

Thus all tests are launched. □ 



